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Abstract 

A comparison is made between conventional Michaelis-Menten kinetics and two 
models that take into account the duration of the conformational changes that 
take place at the molecular level during the catalytic cycle of a monomer. The 
models consider the time that elapses from the moment an enzyme-substrate 
complex forms until the moment a product molecule is released, as well as the 
recovery time needed to reset the conformational change that took place. In the 
first model the dynamics is described by a set of delayed differential equations, 
instead of the ordinary differential equations associated to Michaelis-Menten 
kinetics. In the second model the delay, the discretization inherent to enzyme 
reactions and the stochastic binding of substrates to enzimes at the molecular 
level is considered. All three models agree at equilibrium, as expected; however, 
out-of-equilibrium dynamics can differ substantially. In particular, both delayed 
models show oscillations at low values of the Michaelis constant which are not 
reproduced by the Michaelis-Menten model. Additionally, in certain cases, the 
dynamics shown by the continuous delayed model differs from the dynamics of the 
discrete delayed model when some reactant become scarce. 
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1 Introduction. 



For over a century, the conventional way to repres ent many enzyme react ions 
has b een based on the Michaelis-Menten kinetics ( iMichaelis and Mentenl . 
19131 ) (hereafter MM model), in which substrate S binds reversibly with 
enzyme E to form an enzyme-substrate complex ES, which is later 
transformed into product P and free enzyme E. Derivation of the 



Michaelis-Menten equation is based on the following scheme (ISegell . Il975l ): 



E + S ^ 



^ ES ^ E + P 



However, there is strong evidence suggesting that enzime activity depe nds on 



conformat i onal c hanges in the enzyme structure (IBennet t et al.l. Il978 



Herschlad. 


19881; 


Hakansson et al.. 


I997I; 


Aemon, I2OOO.: .Eisenmesser et al.. 


2OO2I; 


Benkovic and Hammes-Schifferl, 


2003r). It is therefore reasonable to 



assume that a conformational change takes place from the moment an ES 
complex forms to the moment a product P is released; such change requires 
a certain time for its completion. This conformational change must then be 
reset before the enzyme is ready to bind substrate again, a nd so a re c overy 



time must elapse once P is released. Single - enzyme studie s ( iLu et al.l , Il998 



Xie and Lu, 1999; Eisenmesser et al., 2005: Walter , 2006) and nuclear 



magnetic resonance spectrosc opy studies (lEisenmesser et al.l , 12005 : 
Huang and Monteliond , |2005|) support these n otions. A scherne tha t 
considers enzyme regeneration is presented by English et al.l (120061 ), 



E + S ^ 



ES 



fc2 



E" + P, 



E°^E 



where E° represents recovering enzyme. Conventional MM kinetics do not 
account in an explicit way for the times required by these conformational 
changes. 

The aim of this paper is to present two models of enzyme catalysis, one 
continuous and one discrete, that take into account the times needed by 
enzymes to process ES complexes and to reset to its original state once 
products has been released; these models can be considered as an extension 
of conventional MM kinetics which includes delays. The main motivation 
behind the continuous delayed model is to establish the limitations of the 
MM description in non-equilibrium conditions. In order to examine the 
limitations associated to the continuous nature of the first model, we 
consider a discrete delayed model in which the catalytic cicle of each enzyme 
in the system is represented by a re cursive function, hereafter referred to as 
the enzyme map. Previous work by IStange et al.l (119981 , Il999l ) present an 
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analogous description of enzyme action in which each protein molecule is 
represented by a clock-like automaton that binds substrate S and transforms 
it into product P. The time needed by a single automaton to complete this 
transformation and return to its original waiting state represents the enzyme 
turnover time r. The conformational changes of the enzyme (and the 
different stages of the catalytic cycle) are represented by motion along an 
internal 'phase' coordinate. Under the aproppriate conditions this model 
exhibits complex behavior such as oscillations and cluster formation; thus 
another motivation is to explore the emergence of such behavior in our 
model. Unlike other microscop i c models based on stochastic s i mulat ion 



flHaseltine and Rawhnd. l2002l : iKierzekl . l2002l : iRao and Arkinl . 12003 



Puchalka and Kierzekl . 120041 ). our discrete model is completely deterministic, 
and each enzyme is represented by a map instead of a clock-like unit. 



In Section 2 the continuous delayed model is presented by means of a simple 
reaction, and a comparison is made between numerical integration provided 
by MM kinetics and those produced by the delayed model. Section 3 
presents the equivalent discrete model; the results of the discrete simulation 
are compared with results obtained through the continuous delayed model. 
In the conclusions presented in Section 4 we highlight the respective merits 
and limitations of each model. 



2 MM vs continuous delayed model 



To illustrate the effect of delays on the dynamics of enzyme catalysis we 
consider the simple reaction A ^ B, where conversion between two 

substrates A and B is catalyzed by two unidirectional monomeric enzymes a 
and [3. 

As mentioned in the introduction, each conformational change taking place 
along the catalytic cycle will require a certain time for its completion. We 
then distinguish three stages in the cycle: 1) the enzyme is free, waiting to 
bind substrate; 2) the enzyme has bound substrate forming an ES complex; 
3) product P has been released and the enzyme is resetting to its initial 
conformation. The processing time Tp required by stage 2 is a fraction c of 
the turnover time r, while the recovery time Tr associeted to stage 3 is 
(1 — c)r. The turnover time is the time that elapses from the moment an ES 
complex forms until the moment the enzyme is ready to bind substrate 
again; such time is given by 



maxl-^e ^ 
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where Vmax IS the maximum velocity of the enzyme in //mol/ (min-mg of 
protein) and /Xg its molecular weight in Daltons. 

When processing and recovery times are considered the reaction is described 
by a set of four coupled delay differential equations (DDEs): 



^ = ^ {9,mt - ^JW,At - rj - 9AA{t)]M^,{t)} + [AUt)] (2) 
^ = y^ {9.[Mt - rJ]M^,{t - rJ - g,[B{t)]M^,{t)} (3) 
X W-\mt-Tj]M^,{t-Tj - [A{t)]M^,m (4) 



dt 

^ = X 10-% {[B{t - T^)]M^,{t - T,) - [B{t)]M^,{t))} (5) 



where [A] , [B] are substrate concentrations in //M; M^^ , M^^ represent the 
amount of free enzyme in milligrams and yUo,, are the molecular weights of 
enzymes a, j3 in Daltons. The delays and are their respective turnover 
times, whereas r^^ = c^r^, r^^ = c^r^ are their respective processing times. 
is the volume in which the reacion takes place, in liters. The rates and 
are given by g^ = Vmaxi/Ki {i = a,f3), where Ki is the Michaclis constant in 
fiM. [Aext{t)] represents the rate at which substrate A is added to the 
system. Equations for M^2(t)y M^al^); ^/32(^) ^ozi'^) be easily 
obtained from [A(t)], [B(t)], M^^it), M^^{t). For instance, the term 
jjLa X g^[k{t)\M^^{t) is the rate at which M^^{t) appears, while 
^jboc X 10'"^ft«[A(i - T„p)]M,j(t - T„p) is the rate of dissapearance of M^^{t). 

Equations (2-5) reduce to conventional MM equations 



^ 1 (y^^ [^W] M, - Kn„, J^^JL_M I + [K,m (6) 

dt Vr\ ^Kf} + [B{t)] " [A(t)] "J^LextJV^V^ 

rf[B] _ d[k] 
~dr~ "df' 



(7) 



when r = 0. In this case, the amounts of free enzymes iV4^ and M^^ are given 

by 



M 

where and M„ are the total amounts of enzymes a and (3 in milligrams. 
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respectively. Notice that unlike the delayed model (Eqs. 2-5), the MM 
equations (Eqs. 6-7) implicitly assume that substrates are instantaneously 
transformed into products and that the free enzyme amounts are 
instantaneously determined by substrate concentrations. On the other hand, 
in the delayed model the rate at which substrate A appears is proportional 
to the rate at which substrate B was bound by enzyme /3 at time t — r^^. 
Similarly, the rate at which free enzyme (e.g. Afg^(t)) appears at time t 
equals the rate at which free enzyme dissappeared a turnover time r earlier 
(e.g. M,,(t-r)). 

When stationary conditions are reached in the continuous delayed model, the 
amounts of enzyme (in milligrams) in each stage are 



M 

M , = =^ (10) 

" 1 + X 10-3^JAK ^ ' 

u„ X [AIM r , , 

M=— ^ 11 

1 + /i. X 10-3^JA]r 

/X. X 10-\[A]M^r_(l - cj 

1 + X 10-3^JAK ^ ^ 

where the subscripts 1, 2, 3 correspond to free enzyme, occupied enzyme and 

recovering enzyme, respectively, and [A] is the stationary value of 
concentration [A(t)]. The same expressions can be used to determine M^^, 
and M^g by replacing [A] for [B] and the parameters M^, r^, /x^, and 
for those corresponding to enzyme /3. 

In the continuous delayed model the stationary substrate concentration [A] 
corresponds to one of the solutions of a third degree polynomial when 
considering Eqs. (2-5) in stationary conditions. A similar equation must be 
solved to find the stationary concentration [B] . From the solutions of these 
polynomials, the stationary solution is the one satisfaying the restrictions 
< [A] < [So], < [B] < [So], where [Sq] is the total substrate concentration 
in the system. Note that at any time [So] = [A](t) + [B](t) + [a2](^) + [f^2\{t), 
where [02] (t) and [/32] {t) are the concentrations of processing enzymes. 

The corresponding stationary values of the substrate concentrations [A]^j^^, 
[-^]MAf ^'^'^ MM model (Eqs. (6) and (7)) are given by the solutions of two 
second-degree polynomials that satisfy < [A] < [Sq], < [B] < [Sq]. Note 
that at any time [Sq] = [A](t) -t- [B](t). Since the MM description does not 
consider bound substrate, the continuous delayed model and the MM model 
agree in the stationary state when [Aj^^y^ is compared to [A] -|- [0(2]; a similar 
relation applies to substrate B. 

In order to estabhsh a comparison between the delayed model and MM 
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kinetics, numerical integration of Eqs. (2-5) and (6-7) was performed. Initial 
substrate concentrations are set to zero and all enzymes are initialized in 
their idle state. Later, substrate A is quickly added until a concentration [Sq] 
is reached according to [Aextit)] = [Aext,max] exp(-(t - t^Y/a), where 
[Aext,max] IS the maximum rate of substrate addition reached at time t = to] 
the parameter a allows to control the time span during which significant 
amounts of substrate are added. Hereafter we adopt a = O.lr^, and to = 2Ta. 
Unless otherwise stated, the results shown correspond to the following 
parameter values: the reaction volume K - is set to 10.3 attoliters, 



corresponding to a T. brucei glycosome (iNavid and Ortoleval . l2004j ) . For the 
sake of simplicity the reaction is made symmetric, i.e. 
Vmax = Vmax^ = Vmaxp = 500 /xm/(min-mg of protein); the amount of 
enzyme is chosen to be = = 1.66 attograms, a figure equivalent to 
200 molecules of an enzyme with a molecular weight of 5000 Daltons. This 
corresponds to an enzyme concentration [a] = [f3] = 32.24 /iM. For these 
parameter values = = 0.024 sec. In the following example 
Km = Ka = Kp] for the DDE case we use c = = = 0.5. 

Figure la shows a comparison between the evolution of [A] in the MM and 
DDE models for Km = 100 and [Sq] = 100 /iM. When the stationary 
regime is attained, the concentration [A] in the DDE case is below that of 
the MM case since in the DDE model a portion of substrate is bound in the 
form of ES complex with concentration [02] • Figure lb shows results when 
Km is reduced to 0.1 /iM. The continuous thin curve correspond to the 
evolution of [A] + [02]- Several new features are now apparent: a) The 
difference between MM and DDE free substrate concentrations widens 
because in the DDE model the portion of substrate in the form of ES 
complex increases as Km is reduced. And b) oscillations are present in the 
DDE case as a result of the history dependence of the system evolution, an 
usual behavior of systems described by delay differential equations 



(iParravand . Il996l ). Note that a longer characteristic time is needed to reach 
equilibrium when Km is reduced; this is due to the fact that a reduction of 
Km increases the binding rate of both substrates while the fraction of 
available free enzymes decreases, but due to the symmetric cyclic nature of 
the reaction, substrate A is replenished at a faster pace. 



Figure 2 presents the same comparison shown in Fig. lb, but with c lowered 
to 0.02. As stated above, the MM model assumes the products are released 
instantaneously after substrate is bound by the enzime. Therefore, as 
expected, for = ct <^ t substrate concentrations in the DDE model 
oscillate around those in the MM model. 



Figure 3 shows a comparison between the MM and DDE models of the 
fraction of free enzyme fai{t) and the amount of substrate A normalized to 
the total subtrate concentration [Sq] for Km = 0.1 /iM, [Sq] = 100 /iM and 
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c = 0.5. As substrate A is being added to the system the fraction of free 
enzyme fai{t) in the DDE model shows a delay with respect to the MM 
model. This delay is the result of an important difference between MM 
kinetics and the delayed model: in the former, the amount of free enzyme is 
determined instantaneously by the substrate concentrations through Eqs. 8 
and 9; this differs from the time-delayed evolution described by Eqs. 4 and 5. 



3 Continuous delayed model vs a discrete model 



Hess and Mikhailovi (Il994l ) proposed a discrete enzyme model in which each 
enzyme molecule is represented by a clock-like automaton with a phase that 
starts advancing once a substrate molecule binds to it, initiating a catalytic 
cycle whose duration represents the turnover time r. At a certain moment 
along this cycle the automaton releases a product, returning to its initial 
phase after a recovery time has elapsed. Enzyme reactions are thus 
represented in terms of a automata network, in which each automaton is 
coupled to the others through binding and release of intermediate products. 
Instead of a clock-like automaton, the discrete model proposed here assumes 
that the phase x of each single enzyme along the catalytic cycle is 
determined by the recursive function f{x) (hereafter the enzyme map) 
shown in Fig. 4. The map is a deterministic procedure to evolve the phase of 
the enzimes from its phase state Xt at discrete time t to Xt+i at discrete time 
t + 1; that is Xt+i = f{xt). 

The map has a chaotic region defined in < a; < 1, which corresponds to the 
idle state of the enzyme (e.g. ai for enzyme a). Iterations start in this region 
from an initial random phase value Xq] the average number of iterations that 
transcur in the chaotic region is inversely proportional to p, the width 
between the two intersections of the map at / = 1 (p = 0.25 in Fig. 4). 
When an iterate Xt in the chaotic region falls into the range 
{1 — p)/2 < Xt < (l+p)/2 the enzyme is assumed to bind a substrate 
molecule: the next iterate exits the chaotic region and enters the laminar 
region. The laminar region is divided in two sections: the section from x = 1 
to a; = 1 + c6 corresponds to the processing of the substrate by the enzyme 
(i.e. 0:2) whereas the section from x = l + cbtox = l + b corresponds to the 
recovery process of the enzyme (i.e. 0:3). Once Xt > I + cb a. product molecule 
is released, as indicated in Fig. 4. The region of the map between x = 1 + b 
and X = 1 + b + a allows to reinject in a single step the iterate into the 
chaotic region, where the enzyme is ready to bind a substrate again. The 
enzyme map can also be used to model more complicated structures and 
functions such as oligomers and competitive/non-competitive inhibition; 
these are properties that will be presented elsewhere. 
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For a fixed value of a (e.g. a = 1/10), the two parameters of the map b and c, 
and the variable p can be adjusted to obtain an equivalence with the 
continuous delayed model. Note that the number of steps in the laminar 
region is b/a; then if each time step corresponds to a unit time At, the 
turnover time is t — At b/a. The time At associated to each iteration must 
be small enough to produce a significant number of iterations in the laminar 
region; once a given At is chosen the length of the laminar region is set 
according to b = ra/ At. The parameter c = r^/r is the same one defined for 
the continuous delayed model. Variable p is related to the average time ri 
that an enzyme remains in its idle state and is assumed proportional to the 
number Ug of substrate molecules in the simulation volume (i.e. 

Ti = At/p = l/{kins), where ki is the E + S — > ES enzyme rate constant). 
The number of susbstrate molecules that corresponds to a given 
concentration [A] (in /xM) is ua = K-[A]10~^A^^, where Na is Avogadro's 
number. When [A] = Km the velocity of the enzyme is halved; in this 
condition the time spent in the chaotic region must be equal to the time 
spent in the laminar region; that is ri = r. Therefore, p — Atki Ug with 
ki = VmaxfJ'e/ (Q^VrK M^^^^ Na) ■ Finally, the number of enzymes that 
corresponds to a given mass of enzyme rrie in mg is rie = rrieNA/ (/XglO^). 

By using an array of such maps it is possible to consider the effect of the 
discrete nature of reactants on the dynamics of the system and analize the 
conditions under which the continuous delayed model is no longer 
appropriate to describe very small systems. 

Figure 5 shows integration results for the delayed continuous model and its 
equivalent discrete counterpart. In the DDE simulation we use the same 
parameter values as in Figs. 1 and 2, but with Km = K^ = Kf^ = 0.1612 fiM. 
In the system volume of 10.3 attoliters this concentration correspond to 1 
substrate molecule; on the other hand the total substrate concentration [So] 
= 100 fjiM corresponds to 620 substrate molecules. As before, there are 200 
enzyme molecules of each type with a turnover time r = 0.024 s and a 
molecular weight of 5000 Da. Therefore, in the discrete simulation we must 
follow the phase evolution of 400 maps, 200 for each type of enzime. If we 
adopt At — 1 /IS and a — 1/10, the transit in the laminar region occurs in 
24000 iterates (i.e. b — 2400). As the reaction progresses the number of 
available free subtrate molecules {ua and ub) changes, and so the value of p. 
For the maps representing enzime a we have 

Pa = nAX AtVmaxl^e/{QOVrKalO-^NA) = ( At/r) (n^/na) ~ 71^/24000, where 
Ua is the number of substrate molecules corresponding to concentration Ka. 
A similar relation holds for p^. When the n^/ria ratio is 1 the number of 
iterations spent in the chaotic region equals the number of iterations spent in 
the laminar region, so the enzyme will operate at half its maximum velocity. 

As can be seen in Fig. 5, agreement between both models is excellent, but 
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there is one caveat: care must be taken so that the time step At is small 
enough to guarantee that always p < 1; otherwise the results provided by the 
discrete model might disagree from the results of the continuous delayed 
model. In other words, At must be less than TUa/ns, where ns is the 
number of substrate molecules corresponding to the total substrate 
concentration [So]. 

Unlike the discrete model, the DDE approach is expected to fail when the 
number of substrate molecules in the reaction volume is small. As an 
example of this limitation, Fig. 6 shows the evolution of free substrate A in 
both models when [Sq] is reduced to 32.24 fiM, a concentration equivalent to 
200 substrate molecules, half the number of enzyme molecules. The evolution 
of [A] shows how the results provided by the discrete model tend to follow 
those corresponding to the delayed continuous model; however fluctuations 
in [A] can be seen as substrate molecules are bound and released by the 
enzymes. Moreover, synchronization between the delayed continuous model 
and the discrete model is lost as the reaction progresses. This is the type of 
situation that can be expected in conditions of substrate starvation or when 
the reaction volume is reduced while maintaining [So] and the amount of 
enzyme. 

As previoulsy mentioned, in the discrete model it is possible to follow the 
phase evolution of every single enzyme in the system. We show this by 
examining the same simple reaction used so far but with parameter values 
for which sustained oscillations occur; they are K = 5.15 attoliters, 
Vraaxo, = 500 /im/ (mlu-mg) , V„,ax^ = 1000 /im/(min-mg), c = 0.5, = 0.02, 
Ka = 1 A*M, Kp = 0.5 ijM, and [So] = 150 /jM (equivalent to 465 substrate 
molecules) . The number and molecular weight of the enzymes remains the 
same as in previous cases. Figure 7 shows the phase histograms for enzymes 
a for two time moments separated by Ta/4:, ilhistrating how the phases of 
enzimes a tend to be synchronized and how the cluster phase advances 
through the catalytic cycle. In the first histogram there are some enzymes 
waiting to bind substrate, while in the second histogram these enzymes have 
already bound a substrate. 



4 Conclusions 

We have included processing and recovery times in enzyme kinetic models in 
order to account in a simple way for the conformational changes that a 

monomer undergoes during the catalytic cycle and its effects on the resulting 
dynamics. In certain cases, the dynamics exhibited by the continuous 
delayed model differs drastically when compared to that of MM kinetics. 
Examples found in the literature indicate that oscillations in enzyme 



9 



reactions require some strong form of non-linearity: inhibition-activation 
mechanisms, presence of feedback loops, interaction be tween the elem ents of 
a network, or cooperativity between enzyme subunits (ISerkovl. Il968l: 
Goldbetterl. Il996l: iLaub and Loomil llQQSl : iNavid and Ortoleval . 12004 ; 
Chickarmane et al.l . 120071 ). However, the emergence of oscillations in enzyme 
reactions can also be promoted by the delayed dynamics. 

Unlike the MM model, the delayed model accounts for that portion of 
substrate that is present in the system as ES complex; thus stationary 
substrate concentration values agree when the sum of free and bound 
substrate concentrations in the continuous delayed model are compared with 
free substrate concentrations in the MM model. Another important difference 
between MM kinetics and our model is that the fraction of free enzyme in the 
continuous delayed model follows a time-delayed evolution, while in the MM 
case such fraction is instantaneously determined by substrate concentrations. 

The proposed discrete model agrees with the continuous delayed model when 
large enough number of substrate molecules are considered. When this 
number is small the discrete model shows the expected fluctuations that are 
associated to the discrete nature of reactions at microscopic scales. These 
effects can be important when considering reactions that take place in very 
small volumes such as those of organelles. 

The discrete model offers the possibility of following the evolution of each 
single enzyme, and therefore allows for the study of situations in which 
self-organization phenomena arise, such as synchronization and clustering. 
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Fig. 1. a) [A] evolution for MM and DDE models when 

Vmax = Vmaxc, = Vmaxp = 500 /im/(min-mg) , Km = Ka = Kf3 = 100 //M, 
[So] = 100 /^M, rria = mp = 1.66 attograms, c = Ca = = 0.5 and Vr = 10.3 
attoliters. b) Same situation with Km — K^ = Kp = 0.1 /xM (see text). 

Fig. 2. Same as in Fig. lb but with c — Ca — cp — 0.02. 

Fig. 3. Fraction of free enzyme a and normalized concentration [A] for MM 
and DDE models for the parameter values used in Fig. lb. 

Fig. 4. The enzyme map f{x) (see text). 

Fig. 5. Evolution of the amount of free substrate A for DDE model and the 
discrete model. Left hand y-axis shows concentration values while right hand 
y-axis shows corresponding number of susbstrates. Paramater values are 
Vmax = Vmaxa = Vmaxp = 500 /xm/ (miu-mg) , Km = Ka = Kp = 0.1612 /xM, 
[So] = 100 /iM, rria — mp — 1.66 attograms, c — Ca — cp — 0.5, /S.t — 1 
and Vr — 10.3 attohters. 

Fig. 6. Same as Fig. 5 but with [Sq] = 32.24 /xM. 

Fig. 7. Phase histograms for enzyme a. Parameter values are Vmaxa = 500 
/im(min-mg), Vmaxf^ = 1000 /im/(min-mg), = 0.5, Cp = 0.02, Ka ^ 1 /xM, 
Kp = 0.5 /xM, and [So] = 150 /xM and Vr = 5.15 attoliters. 
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